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Atmospheric  wave  data  from  a  multisensor  array  have  been  analyzed. 

It  was  found  that  Saturn  Apollo  launches  generate  acoustic  gravity  waves  which 
propagate  in  the  upper  atmosphere  and  induce  phase-path  fluctuations. 

The  dispersion  study  of  a  TID  event  has  been  used  to  illustrate  the 
possibility  of  deriving  neutral  atmospheric  structure  from  such  type  of 
studies . 

The  use  of  a  parametric  model  to  describe  the  power  spectra  of  the 
phase-path  fluctuations  was  investigated.  The  results  indicate  that  the 
background  activity  in  the  ICPH  to  13  CPH  frequency  range  follows  an  excita¬ 
tion  mechanism  having  12-  and  24-hour  components. 

Phase-path  records  corresponding  to  high  energy  events  were 
studied.  It  was  found  that  identification  of  signals  coming  from  the  events 
is  very  difficult  due  to  the  large  increase  of  the  background  activity  at 
the  expected  signal  arrival  time. 
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1 .  INTRODUCTION 

The  objective  of  this  project  was  to  analyze  atmospheric  wave  data 
collected  with  a  multisensor  array  during  1969-1970.  The  multisensor  array  con 
sisted  of  the  superposition  of  long  period  microbarographs,  magnetometers  and 
phase-path  (doppler)  sounders,  and  was  situated  in  the  New  York-New  Jersey 
region.  The  analyses  have  been  directed,  in  particular,  to  those  events, 

natural  or  man-made,  from  which  atmospheric  structure  information  could  be 
obtained. 

Data  recorded  during  Saturn-Apollo  launchings  have  been  extensively 
analyzed  to  see  if  any  evidence  of  coupling  between  the  upper  and  lower  part 
of  the  atmosphere  could  be  found.  A  study  of  other  rocket  launchings  was  car¬ 
ried  out  to  determine  if  these  rockets  also  produced  doppler  signatures  similar 
to  those  observed  during  Apollo  launchings. 

Some  Traveling  Ionospheric  Disturbances  (TIDs)  were  observed  during 
the  operation  of  the  phase-path  sounder  array.  The  dispersion  study  of  the  TID 
event  of  23  November  1969  has  been  completed  and  is  illustrated  here. 

A  statistical  model  of  the  ionospheric  background  motions  has  been 
constructed.  This  work  was  done  in  cooperation  with  Dr.  Melvin  Hinich  of  the 
Department  of  Statistics,  Carnegie-Mellon  University,  Pittsburgh.  The  model 
describes  the  power  spectrum  of  the  phase-path  fluctuations  as  a  function  of 
frequency  and  time.  The  parameters  of  the  model  have  been  determined  by  multi - 
regression  analyses  of  hundreds  of  power  spectra.  The  effect  of  sampling  on  the 
model  has  also  been  studied. 


tmeleoyne 

ISOTOPES 


A  simplified  theoretical  study  of  the  ionospheric  response  to  interna! 
gravity  waves  was  carried  out.  r„e  approach  to  this  study  has  been  from  the 
point  of  view  of  the  doppler  observations. 

Finally,  the  records  corresponding  to  the  high  energy  events  of  the 
summer  of  1970  have  been  studied  to  look  for  possible  signatures. 
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2.  ROCKET-GENERATED  SIGNATURES  ON  PHASE-PATH  RECORDS 

Introduction 

Phase-path  (doppler)  variations  induced  by  very  low  frequency 

(0,004167  Hz  -  0.0167  Hz)  acoustic-gravity  waves  generated  by  large  rockets 

were  observed  with  an  array  of  phase-path  sounders  at  the  times  of  the  Saturn- 

Apollo  12  and  13  launches.  Array  processing  of  these  signals  showed  that  the 

average  horizontal  phase  velocities  for  these  events  were  of  the  order  of  720 
.  1 

m  sec  and  the  arrival  azimuth  about  175°  (Tolstoy  et  al . ,  1970;  Montes  and 
Posmentier,  1971). 

Since  one  or  two  phase-sounders  and  some  microbarographs  were  also 
in  operation  during  the  launches  of  Saturn-Apollo  8,  9,  10,  and  11,  it  was 
decided  to  review  the  phase-path,  magnetometer  and  microbarograph  records  to 
investigate  if  all  Apollo  launches  have  generated  signals  of  this  kind  and  if 
the  signals  were  present  in  the  microbarograph  and  magnetometer  records.  In 
addition,  we  have  also  searched  the  records  at  times  corresponding  to  other 
rocket  launchings. 

Data  Analyses  and  Results 

The  data  used  in  this  investigation  were  obtained  in  the  New  York- 
New  Jersey  area  with  a  low-frequency  multisensor  array  consisting  of  micro- 
barographs,  magnetometers  and  phase-path  sounders  (Fig.  1).  The  instrumenta¬ 
tion  has  been  described  in  a  previous  report  (Montes  et  al . ,  1970).  The  over¬ 
all  system  operated  with  an  approximately  flat  response  in  the  frequency  range 
0.05  Hz  to  0.0002778  Hz. 

The  phase-path  data  were  recorded  using  phase-path  (doppler)  sounders 
essentially  of  the  type  described  by  Davies  (1962)  but  with  slight  modifications 
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FIGURE  1.  Array  gcoaetry 
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to  permit  automatic  digitization.  The  digitization  process  does  not  permit 
distinguishing  between  the  ordinary  and  extraordinary  modes  of  propagation  of 
the  sounding  radio  waves.  However,  it  has  been  observed  that,  in  general, 
the  ordinary  mode  is  the  dominant  one  and,  therefore,  we  have  assume!  that 
our  observations  are  a  measure  of  this  mode.  The  microbarograph,  magnetometer 
and  phase-path  data  were  filtered  prior  to  digitization,  with  low  pass  electric 
filters  having  the  appropriate  time  constant  to  minimize  aliasing  effects. 

Sounding  frequencies  of  4.824  MHz  and  6.030  MHz  were  used  in  the 
phase  sounders.  At  these  sounding  frequencies  and  at  times  corresponding  to 
the  arrivals  of  the  rocket-generated  signals  the  radio  ref  lection  heights 
were  between  160  km  and  200  km  for  4.824  MHz  and  between  210  km  and  240  km 
for  6.030  MHz.  These  radio  reflection  heights  were  determined  from  true 
height  analyses  of  ionograms  of  the  Hanover  Station  of  Dartmouth  College. 

During  the  Apollo  8  and  9  events  only  one  or  two  phase  sounders 
were  in  operation  and  phase  velocity  estimation  for  these  events  was  not 
possible.  At  the  time  the  signals  generated  by  the  Saturn-Apollo  10  launch¬ 
ing  were  expected  over  the  array,  the  phase-path  and  magnetometer  records  ap¬ 
pear  very  disturbed.  This  appearance  of  the  records  was  due  to  a  Sudden  Iono¬ 
spheric  Disturbance  (SID)  which  started  at  1710  UT  and  lasted  for  about  two 
hours.  The  noise  level  on  the  records  was  so  high  that  they  were  unsuitable 
for  any  type  of  analysis.  If  any  signatures  induced  by  acoustic-gravity  waves 

were  present,  they  were  probably  masked  by  the  noise. 

Doppler  records  for  the  Apollo  11  event  showed  a  high  level  of  noise 
at  some  of  the  array  stations.  This  fact  prevented  the  use  of  an  array  pro¬ 
cessor  on  the  records.  Nevertheless,  a  careful  examination  of  the  Westwood 
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station  records  which  were  the  least  noisy  revealed  that  during  the  time  inter  - 
val  1430  to  1533  UT,  doppler  fluctuations  showing  charact. eristics  similar  to 
those  observed  during  other  launches  can  be  distinguished  on  the  traces  (Fig.  *-)• 
Microbarograph  array  data  Cor  several  of  the  events  have  been  exten¬ 
sively  analyzed  using  frequency  and  velocity  filtering  as  well  as  other  types 
of  filters  and  beamforming  techniques.  A  filter  that  has  proved  to  be  very 
effective  in  other  circumstances,  which  consists  of  cross-correlating  each 
channel  with  the  doppler  trace  containing  the  signal, (same  principle  as  the 
matched  filter)  was  applied  to  the  microbarograph  data.  The  output  of  this 
filter  was  then  beamforined.  This  careful  and  time  consuming  search,  however, 
has  not  been  able  to  detect  any  signals  or  changes  in  character  on  the  pres¬ 
sure  records  corresponding  to  those  found  on  the  doppler  records. 

Magnetometer  recordings  have  been  similarly  analyzed,  although  not 
as  extensively  as  in  the  case  of  the  pressure  fluctuation  data.  The  results 
of  the  analyses  indicated  that  there  were  no  signatures  in  the  magnetometer 
records  at  the  times  that  they  are  found  on  the  phase-path  records. 

Figure  2  shows  the  doppler  records  corresponding  to  the  Westwood 
station.  The  records  have  been  drawn  in  such  a  way  that  the  lift-off  times 
are  lined  up  along  a  vertical  iiine.  The  plotting  scales  are  the  same  for  all 
records.  The  data  have  been  digitally  filtered  with  a  Fourier  filter  having 
a  band  pass  0.0167  Hz  -  0.000333  Hz. 

Visual  inspection  of  the  records  in  Figure  2  shows  the  following: 
a)  At  ab<?ut  56  minutes  after  lift-off  time,  the  character  of  the 

phase-path  fluctuations  changes  from  more  or  less  random-type  small 
amplitude  variations  into  more  coherent  larger  amplitude  wave  trains. 
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FIGURE  2.  Phase-path  records  corresponding  to  Saturn  Apollo  launches. 


"FTTELEDYNE 

ISOTOPES 


8 


b)  The  amplitude  and  frequency  content  of  the  wave  trains  change  with 
time,  suggesting  in  some  cases  a  beating  pattern.  The  maximum  peak- 
to-trough  amplitudes  of  the  wave  trains  are  of  the  order  of  Af  =  0.5 

Hz. 

c)  The  duration  of  this  more  organized  type  of  motion  is  about  60 
minutes  and  the  periods  range  from  about  60  seconds  to  240  seconds. 

d)  Another  feature  that  can  be  recognized  in  the  records  is  a  smaller 
amplitude  (0.2  Hz)  short  wave  train  leading  the  larger  amplitude 
wave  trains.  This  is  indicated  by  a  small  arrow  on  top  of  the 
records  on  Figure  2. 

Power  spectral  analysis  using  ti-e  Blackman-Tukey  algorithm 
(Blackman  $  Tukey,  1959)  was  applied  to  the  unfiltered  data.  The  time  inter¬ 
val  used  in  the  analysis  is  indicated  by  the  arrows  in  Figure  2.  The  results 
of  the  analyses  are  shown  in  Figure  3;  some  common  features  can  be  recognized 
in  the  power  spectral  curves  corresponding  to  the  various  events. 

a)  There  is  a  broad  peak  centered  at  0.0056  Hz  (180  sec)  in  both  the 
4.824  MHz  and  the  6.03  MHz  data. 

b)  The  width  of  the  peak  is  of  about  3  octaves  extending  from  0.00333 
Hz  to  0.0111  Hz. 

c)  The  height  of  the  peak  is  cf  the  same  order  for  all  events.  The 
peaks  of  the  6.03  MHz  data  are  smaller  than  those  of  the  4.824 


MHz  data. 
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Array  processing  of  the  phase-path  data  showed  that  the  arrival  azi¬ 
muth  of  the  signals  measured  at  the  center  of  the  array  (41. TN,  74.12  K) , 
varies  slightly  from  the  average  azimuth  of  l/S*  as  a  function  of  period.  Hie 
largest  change  corresponds  to  the  longest  period  wave  trains  (Fig.  4);  this 
shift  being  eastward  for  the  Apollo  12  and  westward  for  the  Apollo  13  events 
respectively.  The  flight  trajectory  parameters  arc  very  similar  ior  all 
Apollo  launches;  in  particular,  the  altitudes  oi  the  rocket  corresponding  to 
the  various  arrival  azimuths  were  approximately  the  same  for  Apollo  12  and  13 
launchings.  A  major  difference  in  the  trajectory  parameters  is  that  the  sec¬ 
ond  stage  center  engine  was  powered  a  longer  time  (460  sec)  during  Apollo  12 
than  during  Apollo  13  (330  sec);  this  means  that  in  the  first  case  energy  was 


put  at  ionospheric  levels  for  150  seconds  more. 

Assuming  that  the  arrivals  corresponding  to  the  various  launchings, 
for  which  we  do  not  have  array  data,  came  from  the  same  average  directions  as 
those  of  Apollo  12  and  13.  and  using  flight  trajectoty  data  to  locate  the 
source  in  time  and  space,  it  is  possible  to  compute  approximate  .aluvs  of  gm»p 
velocity  as  a  function  ct  period  This  was  dene  by  applying  the  method  ©I 
Ewing  and  Press  (1954)  to  the  record*  shown  in  Pigu.e  2  the  result*  aie  shown 

in  Figure  S. 

Data  corresponding  id  large  rocket  Vomhings  trow  the  let  Unions 


launching  sites  at  Tyuratam  and  Plesetsk  such  Iona  Holniya  It,  Soyu?  i* 
and  others  were  analyzed  to  see  If  some  signals*  could  be  detected  The  re¬ 


sults  of  these  analyses  showed  that  if  the  signatures  were  present,  they  were 
below  the  detection  threshold.  S.m.lar  result*  were  obtained  tar  other  »  hU" * 
launched  from  Cape  Kennedy. 


o 
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Discussion 

The  fact  that  signatures  have  been  identified  on  phase-path  records 
corresponding  to  several  Apollo  launches  confirms  the  hypothesis  that  these 
signatures  are  caused  by  Saturn-Apollo  rockets.  The  average  azimuth  of  arrival 
of  175°  intersects  the  ground  level  projection  of  the  trajectory  at  30.74°N, 
72.76°W.  The  altitude  of  the  rocket  at  this  point  is  160  km.  This  strongly 
suggests  that  the  atmospheric  disturbances  are  generated  at  ionospheric 
heights . 

The  data  presented  here  is  representative  of  winter,  spring  and  sum¬ 
mer  upper  atmospheric  conditions.  The  lack  of  amplitude  and  frequency  content 
variability  from  event  to  event  suggests  that  seasonal  variation  of  lower  atmos¬ 
pheric  conditions  does  not  have  an  important  effect  on  the  propagation  of  the 
acoustic-gravity  waves  that  induce  these  phase-path  signatures.  This  finding 
agrees  with  the  ray  tracing  results  of  Balachandran  and  Donn  (1971)  and  also 
lends  support  to  the  suggestion  made  by  Tolstoy  et  al.  (1970)  that  the  propa¬ 
gation  takes  place  in  the  upper  atmosphere. 

The  failure  to  detect  signals  from  other  rockets  launched  from  Cape 
Kennedy  and  rockets  launched  at  greater  distances  was  probably  due  to  longi¬ 
tudinal  attenuation  of  the  acoustic  waves  or  tc  the  fact  that  the  energy  input 
into  the  atmosphere  at  ionospheric  levels  was  too  small  to  generate  these  long 
period  waves. 

Visual  examination  of  phase-path  sounder  records  for  stations  located 
in  the  Mid-western  United  States  indicates  that  no  signatures  were  present  at 
the  time  of  the  Apollo  12  event  (Davies,  1970),  and  Apollo  15  (Lerfald,  1971). 
This  lack  of  arrivals  could  be  due  to  attenuation  since  the  distance  from  the 
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source  to  these  stations  is  more  than  twice  the  distance  from  the  source  to 
Westwood,  However,  it  is  also  possible  that  this  is  due  to  the  better  response 
of  the  ionosphere  to  neutral  atmosphere  waves  traveling  in  an  almost  North- 
South  path  (Hooke,  1970>  as  is  our  case. 

The  difference  in  arrival  azimuth  for  the  longer  period  signals  dur~ 
ing  Apollo  12  and  13  events  could  be  dye  to  the  longer  duration  of  the  powered 
stage  on  Apollo  12;  since  this  longer  duration  mode  the  rocket  exhaust  plume 
travel  more  towards  the  east  than  it  did  on  the  Apollo  13  event.  Khy  the 
longer  periods  are  the  most  affected  by  this  situation  is  something  we  cannot 
answer  at  this  time.  In  order  to  answer  this  question  one  would  have  to  solve 
the  moving  source  problem,  mis  is  complicated  because  the  component  of  the 
source  velocity  in  the  direction  of  the  array  results  in  appreciable  dopplcr 
shifting  of  the  source  spectrum,  while  the  vertical  component  of  motion  intro¬ 
duces  a  modulation  of  each  mode.  The  latter  effect  is  probably  small  since 

the  change  in  altitude  of  the  rocket  during  the  generation  of  the  waves  Is  not 
too  large. 

Figure  5  shows  that  the  group-velocity  points  are  grouped  into  four 
separate  groups.  The  group  whose  -eio.:ty  and  period  ranges  .ire  450>tl>400  and 
240>T-200  respectively  has  been  identified  by  T:i>toy  ct  al .  (1970)  as  a  mode 
whose  energy  is  trapped  in  the  lewet  ther»Sph«e.  live  other  groups  probably 
belong  to  modes  which  are  trapped  in  the  upper  Saw,d  channel. 

Finally,  it  should  be  pointed  out  that  absence  of  signatures  from 
smaller  rockets  and  the  wavelengths  involved  in  these  observations  strongly 

suggest  that  the  acoustic -gravity  wa/es  are  generated  by  the  rocket  exhaust 
rather  than  by  a  ballistic  mechanism 


3.  DISPERSION  STUDV  OF  THU  TID  EVENT  OF  23  NOVEMBER  1969 


Introduction 

The  earth's  atnosphcre  is  a  v?ry  complex  nediun  and  there  are  severe 
gaps  In  the  knowledge  of  atmospheric  parameters  and  structure,  In  particular 
for  altitudes  above  100  kn. 

In  the  lower  100  la  of  the  atmosphere  the  structure  is  reasonably 
well  known,  and  although  vho  teaperutute  and  winds  \ ary  with  season  and  geo¬ 
graphical  location  it  is  possible,  in  principle,  to  include  these  effects  In 
propagation  studies  (Pierce,  I960;  Balachandrnn,  ID^O).  However,  it  has  been 
observed  that  following  atmospheric  nuclear  explosions,  long  period  (2-5  min) 
acoustic  waves  and  internal  gravity  waves  propagate  at  heights  above  200  km 
(Baker  and  Davie*,  1968;  Obayashi,  1963;  Herron  and  Montes,  1970).  These  ob¬ 
servations  indicate  that  appreciable  wave  energy  in  the  acoustic  and  gravity 
wave  frequency  range  corresponding  to  periods  longer  than  a  minute  can  be 
channeled  at  mesopause  and  theinospherc  altitudes. 

It  is  well  known  that  in  wave  systems  geneiated  by  atmospheric  ex¬ 
plosions  the  distribuMon  of  »sdc  encigy  with  he.ght  depends  on  the  utnos- 
phetic  neutral  temperature  sttuciutc 

Calculations  of  yield  and  height  of  bu.'st  based  on  mode  analysis 
arc  at  the  present  time  marginal  because  of  the  vatiation  of  the  temperature 
structure  Tolstoy  and  Uu  (19*0)  have  shown  that  any  ptediction  of  the  ver¬ 
tical  amplitude  distribution  of  displacement  c(a),  is  prevented  by  the  vari¬ 
ability  of  the  neutral  structure  since  the  displacement  CU)  is  heavily  depen- 
dent  on  the  Biunt-V3ts3U  N(i)  and  acoustic  cutoff  «&  (c)  frequencies. 

A  study  of  natural  ionospheric  backgeound  mot  ons  made  by  Tolstoy 
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and  Montes  (1971)  has  shown  that  it  is  possible  to  derive,  in  situ,  long  term 
averages  of  the  buoyant  resonance  frequency  (VKisdlS  frequency)  of  the  neutral 
gas.  This  is,  of  course,  far  from  the  precise  determination  of  N(z)  and  wo^z) 
that  one  would  like  to  have. 

Tolstoy  (1970)  has  pointed  out  that  dispersion  studies  of  Traveling 
Ionospheric  Disturbances  (TIDs)  could  provide  some  information  about  the  neu¬ 
tral  structure  at  ionospheric  heights. 

During  the  operation  of  phase-path  sounder  array  by  Teledyne  Isotopes 
several  TIDs  were  observed.  Their  general  properties  have  been  reported  by 
Rao  (1970).  In  order  to  test  Tolstoy's  suggestion,  a  TID  was  selected  and 
subjected  to  dispersion  analysis.  The  results  of  this  study  are  reported  here. 

Data  and  Method  of  Analysis 

The  TID  event  selected  for  this  study  is  shown  in  Figure  6.  At  the 
time  of  the  event  there  were  five  doppler  sounders  m  operation.  The  four 
sounders  forming  the  array  operated  at  a  sounding  frequency  of  4.824  MHz; 
the  fifth  sounder  operated  at  Westwood  using  a  sounding  frequency  of  2.412  MHz. 
The  magnetograms  corresponding  to  the  time  interval  of  the  TID  event  showed 
very  little  activity.  This  fact  can  be  taken  as  an  indication  that  the  dis¬ 
turbance  observed  on  the  phase-path  records  had  a  hydrodynamic  origin,  i.e., 
it  was  induced  by  an  atmospheric  grayity  wave. 

The  starting  time  of  the  event  was  set  at  2000  UT  on  November  23, 
I9b9.  However,  it  is  more  difficult  to  estimate  the  ending  time  of  the  TID. 

If  too  long  a  record  is  included  in  the  analysis,  the  doppler  fluctuations 
which  are  not  connected  with  the  TID  will  make  the  coherency  drop  and  will  in¬ 
troduce  errors  in  the  azimuth  and  phase  velocity  estimates.  On  the  other  hand, 
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if  the  length  of  the  record  chosen  is  too  short,  one  could  be  taking  out  of 
the  analysis  important  contributions  to  the  Til)  and  the  frequency  resolution 
of  the  analysis  will  suffer.  With  the  help  of  cross  spectral  analyses  the 
ending  time  for  the  event  was  chosen  to  be  at  0130  UT  on  24  November  I960 

The  dispersion  analysi-  of  the  phase-path  records  at  the  time  Inter¬ 
val  2000  UT  -  0130  UT  was  carried  out  using  the  cophase  analysis  described  by 
Posmentier  and  Hermann  (1971).  The  results  of  the  analysis  have  been  plotted 

as  points  on  the  w-k  plane  (Fig.  7).  Two  conclusions  can  be  drawn  from  this 
figure: 

a)  the  points  fall  on  the  phase  velocity  range  120  to  240  ■  see”1 

b)  they  probably  belong  to  various  modes. 

The  relative  phases  as  function  of  frequency  between  the  2.412  MHz 
data  and  the  4.824  MHz  at  Westwood  were  found  from  the  cross  spectral  analysis 
of  these  two  records.  The  reflection  heights  corresponding  to  the  two  sound¬ 
ing  frequencies  were  obtained  from  the  true  height  analysis  of  the  Hanover 
station  of  Dartmouth  College.*  The  average  reflection  heights  were  found  to 
be  248  km  for  the  4.824  MHz  and  233  km  for  2.412  MHz  respectively. 

Discussion 

The  points  on  Figure  7  can  not  be  joined  arbitrarily.  A  propaga¬ 
tion  model  must  be  constructed  assuming  reasonable  values  for  the  parameters 
until  a  fit  to  the  data  can  be  obtained. 


* 


Made  available  to  us  through  the  kindness  of  l)rs.  M.  G.  Morgan  and 
C.  Calderon  of  the  Radio  Physics  Laboratory. 
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Tolstoy  and  Uu  (1970)  hive  shown  that  the  region  of  the  atmosphere 
rtete  »0<  N  must  be  taken  into  account  when  studying  the  propagation  of 
atmospheric  waves  is  above  100  km.  Following  these  authors,  we  write  the 
solution  of  the  differential  equation  for  the  pressure  as 


P  *  o'1^2  q  exp  [(ax  -  wt)] 


(1) 


inhere 

o  is  the  neutral  gas  density, 
a  Is  the  horizontal  wave  number, 
u  Is  the  angular  frequency  of  the  disturbance  and 
q  satisfies  the  lie  Inholt:  equation 


q..  ♦  6“q  *  o  ; 


(2) 


the  sub-inde.  :  indicates  d.t (crepitation  with  respect  to  that  variable 
and  3  is  a  function  of  :  and  is  given  by 


3“ 


with 


and 


if 


lfi2L2 
4*  c4 


f 


(t-D 


S.  i* 

c2* 


(lnC2)| 


(4) 
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wnere  C  is  the  speed  of  sounu,  g  the  gravity  acceleration  and  t  the  ratio 

of  specific  heats  is  taken  to  be  1.4. 

The  variation  of  the  speed  of  sound  with  z  above  certain  », 'light 

can  be  approximated  by  assuming  a  law  for  C(z)  of  the  form. 

C2  (z)  =  Kz  +  Co2  , 

where  C0  is  the  sound  velocity  at  z  =  0,  and  K  is  a  parameter  to  be  deter¬ 
mined  once  the  atmospheric  model  is  assumed. 

From  equations  (4) ,  (5)  and  (6)  we  obtain  ^  and  N  in  function  of  the 
new  parameter  K 

„  2  U)  -  1  CgV  -  K2)  <7) 

op 

and 

N2 (z)  *  (K  ♦  0.4g)  .  (8) 

<r 


If  now  the  atmospheric  structure  is  approximated  by  taking  the  origin  of 
coordinates  at  h  «  120  km,  the  parameter  K  can  be  determined  from  equation 
(6)  using  values  for  C0  and  C  (140)  from  the  U.S.  Standard  Atmosphere  tables 

C*  =  3.92xl02  m  sec*1  z  =  0 

2  1  5 

C  =  9.90x10  m  sec'  2  *  1.40x10  m 


This  gives 


K  =  5.92  m  sec 


(9) 


and 


-2  -1 
N  =  2.40x10  rad  sec 

_2 

a)  =  1.46x10  rad  sec 
op 


at  z  =  O 


(10) 
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-2  -1 
N  =  1.18x10  rad  sec 


\  at  z  =  1.40x10  m 

Since  the  observations  presented  here  belong  to  altitudes  of  220  - 


-2  -1 
w  =  0.57x10  rad  sec 
op 


(H) 


250  km,  and  we  are  not  concerned  with  amplitude  calculations  at  lower  levels, 
we  do  not  have  to  worry  about  the  structure  below  120  km  and  this  can  be 


approximated  in  the  form  shown  in  Figure  8. 

Using  equations  (3),  (6),  (7)  and  (8)  with  (9),  values  of  3  were 
calculated  for  the  various  points  of  the  w-k  plot  (Fig.  9). 

It  has  been  pointed  out  in  the  infrasonic  literature  that  when  in¬ 
terpreting  observations  of  TIDs  special  care  should  be  taken  to  see  if  these 
observations  belong  to  a  turning  point  (Tolstoy,  1971).  It  is  well  known  that 
WKB  approximations  are  not  valid  at  turning  points,  and  the  interpretation  of 
the  w-k  data  should  be  done  taking  into  consideration  this  fact.  To  test 
our  data  for  proximity  to  a  turning  point  we  have  used  the  criteria  given  by 
Tolstoy  (1971). 

If  —  =  ~  far  from  a  turning  point,  (12) 

Az 


and 


if 


0)  CO 

—  >  >  — 

8  A<j> 


Az 


near  a  turning  point. 


03) 


where  A<j>  is  the  phase  difference  between  the  observations  made  at  two  differ¬ 
ent  levels  and  Az  is  the  vertical  distance  between  these  levels.  The  results 
of  the  cross-spectral  analysis  between  the  4.824  and  2.412  MHz  data  showed  that 
the  phase  differences  are 
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14  *  51 ,5  :  5*  lor  -  *  0  002221  rad  sec"1 

Aft  *  51.0  i  4*  w  *  0.0031*55  Mil  >cc  * 

Aft  *  52.2  2  ?*  «  *  0  0044*3  rid  sc/1 

the  corresponding  vjlur^  of  »i/(£*/3i)  were  F.O  s  sec  *»  53.5  •  sec  and 
73.2  n  iec‘*  respectively.  Hie  values  of  <*/&  obtained  for  the  same  angular 
l tetpten  i les  were  20.0  n  sec  45  O  n  sec  ^  and  50.37  m  sec  respect nely. 

It  can  be  seen  that  taking  into  account  the  eaperimental  errors,  the  data 
shown  here  satisfies  condition  (12)  indicating  that  our  observations  have  been 
node  far  frm  a.  turning  point . 

A  refined  value  ol  S  can  be  obtained  if  we  set  and  use 

this  adjusted  value  of  i  to  recompute  ,V  and  w.  After  this  was  done  »-k 
curves  using  £  as  a  parameter  were  computed  and  the  results  are  shown  In  Figure 
0.  A  study  of  this  figure  indicates  that  our  earlier  conclusion  about  the 
points  on  the  «-k  plane  was  correct,  i.c  ,  the  points  <o  respond  to  various 
modes 

ft  has  been  i ilust'ated  here  that  **  .a  possible  to  derive  structure 
information  from  the  dispersion  ol  naturally  ©.  iiimug  TIPs.  Continuation  of 
this  type  of  study  would  secute  more  data  about  the  neutral  tenpetature  sttuc- 
tore  at  ionospheric  altitude!. 


4. 


PARAMETRIC  MODEL  OF  IONOSPHERIC  BACKGROUND  MOTIONS 


Current  studies  of  atmospheric  acoustic  gravity  wave  propagation  have 
shown  that  there  is  a  wide  distribution  of  the  energy  of  the  various  modes 
throughout  the  atmosphere,  including  the  ionosphere.  Consequently,  it  was  nat¬ 
ural  to  think  that  observations  of  ionospheric  motions  produced  by  motions  of 
the  neutral  component  of  the  atmosphere  could  be  used  to  detect  acoustic-gravity 
wave  modes  generated  by  nuclear  explosions  and  propagating  at  ionospheric 
heights.  Disturbances  associated  with  nuclear  explosions  in  the  atmosphere 
have  indeed  been  recorded  on  both  ground  level  microbarographs  (Donn  and  Shaw, 
1967;  Tolstoy  arid  Herron,  1970)  and  at  ionospheric  heights  (Obayashi,  1963; 

Baker  and  Davies,  1968;  Herron  and  Montes,  1970). 

A  simple  method  of  monitoring  ionospheric  motions  is  to  measure  the 
variations  of  the  phase-path  of  an  ionospherically  reflected  H.F.  radio  wave 
(Davies,  1962).  In  1969,  with  this  objective  in  mind.  Teledyne  Isotopes  set 
up  a  four  element  phase-path  sounder  array. 

It  is  well  known  that  in  order  to  describe  the  performance  of  an 
array,  the  knowledge  of  the  background  characteristics  of  the  parameter  to  be 
measured  is  essential. 

Power  spectral  analysis  is  often  used  to  investigate  the  character¬ 
istics  of  the  background.  The  power  spectrum  of  the  phase-path  variations  is 
a  function  of  several  parameters.  It  reflects,  in  particular,  the  spectral 
characteristics  of  the  disturbance  sources,  the  manner  in  which  this  mechanism 
is  coupled  into  the  medium,  the  efficiency  of  coupling  to  changes  in  the  elec¬ 
tron  density,  and  the  structure  of  the  medium. 
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Ionospheric  background  motions  are  known  to  vary  over  time  (Georges, 
1967;  Tolstoy  and  Montes,  1971).  This  variation  can  be  appreciable  at  times, 
with  large  changes  occurring  over  a  period  of  several  hours.  Computation  of 
the  long-term  average  spectrum  indicates  how  the  power  is  distributed  as  a 
function  of  frequency  but  is  of  little  use  for  predicting  the  background  mo¬ 
tion  for  any  particular  time  of  day.  The  objective  of  this  work  was  to  make 
a  realistic  and  physically  compatible  parametric  model  for  the  power  spectrum 
of  the  phase-path  fluctuations  in  which  the  variation  of  the  parameters  over 
time  accounts  for  the  nonstationarity  of  the  background  motion. 

Data  and  Method  of  Analysis 

During  the  years  1969  and  1970,  Teledyne  Isotopes  recorded  micro- 
barographic,  geomagnetic  and  phase-path  fluctuation  data  in  the  New  Jersey 
area  as  part  of  a  program  to  study  atmospheric  background  motions.  The  set 
of  data  used  in  this  study  corresponds  to  the  month  of  June  1969  (Fig.  10), 
and  it  was  chosen  because  of  the  author's  familiarity  with  it  from  a  previous 
study.  Since  our  main  concern  is  the  phase -path  (doppler)  sounder  data,  these 
will  be  described  briefly.  The  sounder  operated  at  almost  vertical  incidence 
(separation  between  transmitter  and  receiver  was  less  than  30  km),  with  a 
sounding  frequency  of  4.82  MHz;  at  this  frequency  the  ionosonde  was  effectively 
sampling  the  F-regicn  of  the  ionosphere.  The  data  were  automatically  digitized 
after  being  filtered  with  a  120  second  low  pass  filter,  using  a  sampling  inter¬ 
val  of  60  seconds,  in  this  way  minimizing  aliasing  effects. 

The  digitizing  process  does  not  distinguish  between  the  ordinary  and 
extraordinary  modes  of  propagation  of  the  sounding  radio  wave;  since  it  is 
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FIGURE  10a.  Phase-path  fluctuation  data  used  in  the  parametric 
model. 


29 


V*^,A,vW^ 

VyivjkKV^ 

«^w 

v'~;Wr^VYw^A^i/^^'  ' 


"AwA 


**f  >*>»  v  » 


•#?*«* 


j  V  - 

VW'*/“>vv ArW^Vw'*''"’'^^ 


,’Vfr  OfiOi  c*i>6  i'll*  •*» 

ul 


MQim  lob.  continuation  of  Figure  10a. 
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commonly  observed  that  the  ordinary  mode  overshadows  the  extraordinary  one,  we 
have  assumed  that  our  data  provides  a  measure  of  the  ordinary  mode. 

The  data  sample  representing  648  hours  (27  days)  was  subjected  to 
power  spectral  analysis,  and  the  Log^  °f  the  power  spectral  density  was  com¬ 
puted  for  every  three-hour  interval  using  the  Blackman -Tukey  algorithm 
(Blackman  and  Tukey,  1959).  Special  care  was  taken  to  remove  from  the  analy¬ 
ses  those  time  intervals  corresponding  to  solar  flares  or  magnetic  storms 
since  these  events  are  not  part  of  the  background. 

In  order  to  test  the  effect  of  sampling  on  the  analysis,  2,  3,  4  and 
6-hour  time  intervals  were  used  in  the  computation  of  the  power  spectra  for  a 
short  data  sample  (240  hours).  It  was  found  that,  in  general,  the  length  of 
the  time  interval  did  not  greatly  affect  the  amount  of  variance  that  can  be 
explained  by  the  model.  However,  taking  into  account  aliasing  problems  and 
the  desired  frequency  resolution,  it  was  a  natural  compromise  to  arrive  at  the 
three-hour  sampling  intervals  used  throughout  the  analysis. 

The  frequency  band  covered  by  the  spectral  analysis  was  divided  into 
13  elementary  frequency  bands  of  equal  bandwidth  Af  -  0.0002778  Hz.  These 
elementary  frequency  bands  correspond  to  the  center  frequencies  1.0,  2.0, 

3.0, . 13.0  CPH,  that  is,  periods  of  60,  30,  20, . 4.61  minutes,  respective¬ 

ly.  Typical  power  spectra  are  shown  in  Figure  11. 

The  results  of  the  power  spectral  analysis  indicate  that  the  iono¬ 
spheric  background  motions  constitute  a  nonstationary  process,  i.e.,  the  power 
and  the  spectral  shape  of  the  phase-path  fluctuations  vary  over  time  (Fig.  12). 

The  average  shape  of  the  phase-path  fluctuation  spectrum  suggests 
that  a  linear  approximation  of  the  log  spectrum  could  be  used  as  a  model.  The 


FREQ  (C/H) 


32 


33 


statistics  of  such  type  of  approximation  have  been  given  by  llinich  (1967). 

Let  S|j(<fc)  denote  the  average  power  in  the  band  of  center  frequency 

ffc  computed  from  the  Jth  3-hour  long  record  on  the  Ith  day  (1*1 . . 

j*l,2 . 8;  and  k*l,2 . ,13)  and  let  Sj(f^)  denote  the  true  power,  i.e., 

the  long  term  average  value  of  Sjjff^)  **  S  ( f ^ ) .  It  Is  well  known  that 

* 

the  variance  of  Logi0  S^ffj.)  is  the  same  for  all  frequencies  ffc  (Hinlch 
and  Clay,  1968).  the  sampling  variability  of  LogJ0  s|jCfk)  about  the  true 
value  Log|0Sj(f^)  depends  only  on  the  record  length  and  on  the  bandwidth. 

Suppose  that  we  moke  the  following  linear  model  for  Lo8|0Sj 5 
U»I0  W  ■  “  *  *Jfk  *  FJ 

where  a  is  a  constant  effect,  8j  is  the  siope  of  Log|0S^(fk)  as  a  function 
of  ffc  for  a  fixed  J,  i.e.,  we  are  assuming  that  at  any  given  time  of  day  the 
log  spectrum  of  the  background  is  linear  in  frequency.  Fj  is  the  effect  due 

to  the  lime  variation,  this  linear  model  yields  the  following  expression  for 

* 

Log  jo  **ij^|^*  ^or  i  *  l»  •  •  •  *  #n 

U*I0  *IJtfk>  *  1  •  Vk  •  f)  '  Cl,rk' 

where  tjlf^)  is  a  sero  mean  Gaussian  error  whose  variance  depends  only  on  k. 

It  is  clear  that  for  the  model  to  have  validity,  the  time  variation 
effect  oust  be  modeled  taking  into  account  actual  physical  situations. 

During  the  last  two  decades  considerable  effort  has  been  put  in  under¬ 
standing  the  upper  atmosphere  dynamics.  It  Is  in  general  agreement  among  upper 
atmospheric  workers  (Bolgiano,  1959),  that  although  turbulence  is  a  common 
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occurrence  in  the  atmosphere,  at  least  to  an  altitude  of  100  km,  the  intensity 
of  fully  developed  turbulence  may  be  no  more  than  1  or  2  percent  of  the  total 
energy;  the  bulk  of  the  energy  probably  residing  in  more  organized  type  motions 
such  as  random  gravity  waves,  thermal  winds,  and  tidal  oscillations.  Lindzen 
(1969)  pointed  out  that  the  most  significant  atmospheric  tides  are  those  ex¬ 
cited  by  the  absorption  of  insolation  by  water  vapor  and  ozone.  These  tides 
have  periodicities  of  12  and  24  solar  hours.  In  the  upper  atmosphere  the  pic¬ 
ture  is  more  complicated  because  of  the  presence  of  ioni##|ion. 

Simplified  theory  of  ionospheric  layer  formation  predicts  vertical 
motion  of  the  ionosphere  as  a  function  of  solar  angle  (Rishbeth,  1963).  Thus, 
we  have  assumed  as  a  starting  point  for  modeling  the  time  variation  effect, 
a  solar  control  of  the  ionospheric  motions.  Sinusoidal  functions,  with  argu- 
ment  of  the  form  ut  *  .t,  were  generated.  The  parameter  t  corresponds  to 

the  time  at  the  center  of  each  3-hour  window  of  the  power  spectral  analysis,  and 
T  is  a  multiple  or  submultiple  of  24  hours. 

In  order  to  investigate  the  relationships  between  the  sinusoidal  time 
functions  and  the  Log^  power  spectrum  for  each  frequency  band,  multiple  linear 
regression  analysis  was  used.  The  sinusoidal  functions  entered  the  analysis  as 
independent  variables,  while  the  LogjQ  power  spectrum  was  the  dependent  variable. 
Means,  standard  deviations,  simple  and  multiple  correlation  coefficients,  etc. 
were  calculated. 

Results 

The  multiple  regression  analysis  showed  that  the  stronger  contribu¬ 
tions  to  the  time  variation  of  the  power  spectra  come  from  those  periodic  func¬ 
tions  with  angular  frequencies  a  *  2ir  rad  day"1  and  w=  4  it  rad  day’1.  Figure  13 
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FIGURE  13.  Amplitude  and  phase  for  12-  and  24-hour  components  of  time  variation  as  function  of  frequency. 
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shows  the  distribution  of  amplitudes  and  phases,  as  function  of  spectral  fre¬ 
quency,  for  these  functions.  It  can  be  seen  in  this  figure  that  both  the 
amplitude  and  the  phase  of  the  24 -hour  component  (w  =  2rc  rad  day'1)  are  approxi¬ 
mately  constant  for  frequencies  higher  than  6CPH.  A  large  change  in  the  phase 
between  3  and  5  CPH  is  also  observed.  The  amplitude  of  the  12-hour  component 
is  rather  constant  from  2  to  13  CPH.  The  amplitudes  of  the  other  periodic 
components  are  smaller  by  an  order  of  magnitude. 

The  average  slope  of  Log^  (f^)  is  given  by  the  slope  of  the  inter¬ 
cepts.  Figure  14  which  is  a  plot  of  the  intercepts  as  function  of  frequency 
chows  that  the  slope  changes  at  a  frequency  of  7.25  CPH.  The  values  of  B 
corresponding  to  these  two  branches  are 

8i  *  -0.188  «  7.25  CPH 

Sj  -0.066  fR  >  7.25  CPH  . 

Examination  of  the  multiple  correlation  coefficients  reveals  that  the 
model  can  explain  a  good  part  of  the  time  variation  of  the  power  spectra  for 
the  frequencies 

f.  <  3  CPH 
k 

and 

f.  >  6  CPH: 
k 

i.e.,  between  36%  and  50%  of  the  variance  can  be  explained  at  these  frequencies. 

The  distribution  of  the  residuals  for  each  frequency  band  is  shown  in 
Figure  15  in  the  form  of  smoothed  histograms.  The  normal  distribution  curve  is 
also  plotted  in  each  case. 


FREQ  ( CPH) 

FIGURE  14.  Plot  of  intercepts  as  function  of  frequency. 


Mien  the  data  was  ©token  into  a  couple  of  pieces,  stww  indication  o* 
noimationarity  ***  tound  llnwcer  this  fact  »s  «» i  to©  significant  since  one 
mould  anticipate  tut  lor  longer  tecotJs  ^we  seasonal  variation  is  introduced. 

Conclusion 

Iron  the  results  one  can  conclude  that  ionospheric  background  actions 
in  the  frequency  range  coveted  by  the  analysis  resfwwid  to  an  excitation  ntekhan* 
which  has  a  dm  in  cstqnment  of  -■*  Hour?  (diurnal!  and  a  sail  ter  one  of  12 
hoots  (semidiurnal ). 

Ihc  reason  shy  the  node!  cannot  explain  nwch  of  the  variance  In  the 
ireqoency  range  i  Oil  >  f  *»  ©CHI  is  probably  that  this  range  also  corresponds 
to  the  itutn»vaisSI3  frequency  at  ionospheric  altitudes.  It  is  sell  known  that 
the  VSisliS  frequency  besides  being  the  cutoff  frequency  for  internal  gravity 
cases  is  also  the  local  value  of  the  resonant  buoyancy  frequency,  and  It  should 
be  expected  that  locally  generated  disturbances  wilt  show  up  strongly  in  this 
frequency  range. 

So  particular  patterns  have  been  found  In  the  distribution  o«  the 
residuals  -  I  ron  analysis  ol  saoushed  histograms  computed  from  the  residual* 

•hey  ajipeit  to  be  normally  distributed 

finally  it  mist  be  pointed  oo»  that  a  sis  lar  analysis  earned  cut 
u%mg  data  fra®  a  long  period  mi crobai ©graph  shewed  that  no  wore  than  9*  'l 
the  » at  lance  could  be  explained  by  the  model . 


S.  IONOSPHERIC  RESPONSE  TO  INTERNAL  GRAVITY  WAVES 
RELEVANCE  TO  DOPPLER  OBSERVATIONS 

The  response  of  the  F  region  of  the  ionosphere  to  interna]  gravity 
waves  is  anisotropic  due  t;  the  presence  of  the  geomagnetic  field.  Hooke 
(19/0)  has  shown  that  the  response  of  the  ionosphere  besides  being  aniso¬ 
tropic  is  also  subjected  to  diurnal  and  seasonal  changes.  This  author 
describes  the  response  in  terms  of  electron  density  variations,  which  are 
not  directly  related  to  the  dopplcr  observations.  Since  the  phase-path 
sounders  are  assumed  to  measure  the  vertical  velocity  of  the  electrons,  it 
was  decided  to  make  a  simplified  model  of  the  response  in  terms  of  this 
parameter  to  aid  in  the  interpretation  of  our  observations, 

Bugnolo  (196?)  has  shown  that  the  vertical  velocity  of  the  elec¬ 
trons  when  no  external  electric  fields  are  present  can  be  expressed  in 
terms  ol  the  velocity  of  the  neutral  gas  by  the  uncoupled  equation 


V:c  ~  -  Cos  (a)  Cos  (b)  .Cos  ($^)  .U*Cos~  (b)  ,W 

whete  \zc  is  the  vertical  velocity  ol  the  electrons;  the  parameters  refer 
to  Figure  16  and  are 

a  •  angle  between  magnetic  south  and  iio, 

b  angle  between  the  vertical  axis  of  coordinates  and  Ho, 

•j  -  angle  of  arrival  of  the  neutral  wave  measured  from 
magnetic  north, 

Ho  *  the  main  geomagnetic  field, 

U  hor i rental  particle  velocity  of  the  neutral  gas, 

W  a  vertical  particle  velocity  of  the  neutral  gas 


FIGURE  16.  Geometry  for  ionospheric  response. 


In  the  northern  hemisphere 


a  =  Il-I 

b  =  ir/2*I 


(15) 


where  I  is  the  magnetic  dip  angle. 

The  values  of  U  and  W  can  be  obtained  from  equations  (4.2)  and  (4.8)  of 
Tolstoy  (1963)  in  the  form 


U 

W 


-  ia 

2  /r2 
uj  /C  -a 


2 


(v2~g/C>iY 


(16) 


where  the  sign  in  front  of  y  indicates  upgoing  or  downgoing  waves,  and 
a  =  horizontal  wave  number, 

<i>  =  angular  frequency  of  the  internal  wave, 

C  =  speed  of  sound, 

v  =  1/2H,  H  being  the  scale  height, 


and 


2 

Y 


2 
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2 

a" 
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(17) 


with 


N  =  V'aisala  frequency 


Using  equations  ( i 5 )  and  (16)  the  response  of  the  ionization  can  he  written 
in  the  foxin 


V 

.,  zo  2 

R  k - a  -Cos  I.  S  in  I  .host  n.l->Sin  I 

W 


(IK) 


Where  F  =  U/W. 


43 


It  should  be  pointed  out  that  R  is  a  cample*  t.unbcr  and  a>  such  it  has 
amplitude  and  phase. 

Taking  typical  values  Tor  the  vanous  parameters  and  noting  that 
the  magnetic  dip  angle  at  ifestwoed  is  ?•%  one  obtains  the  responses  shown 
in  Figures  1?  through  21. 

Inspection  of  the  figures  indicates  that  for  horizontal  wave¬ 
lengths  of  400  km  there  is  a  large  phase  shift  between  the  electron  vertical 
velocity  and  the  neutral  gas  vertical  velocity  for  arrivals  having  azimuths 
between  -60*  and  60*  while  the  amplitudes  are  about  a  factor  of  three  smaller 
for  arrivals  coming  from  the  south.  Smaller  wavelengths  are  less  affected 
by  this  anisotropy.  The  amplitude  response  curves  compare  well  with  those 
of  llooke  (1970).  Inclusion  of  the  solar  control  on  the  production  of  ioni¬ 
zation  will  distort  these  response  curves  in  the  sane  way  as  those  shown  by 
llooke  The  curves  shown  hen  would  cot  respond  to  noontime  conditions. 
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FIGURE  18. 

Amplitude  and  phase  response 
for  Xx*400  Xz=275  T*26.5 
dip®72° 
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Amplitude  »ml  phose  response 
for  Xm4O0  U*275  T»26.5 
dip*72* 
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6.  ANALYSES  OF  RECORDS  CORRESPONDING  TO  HIGH-ENERGY  EVENTS 

During  the  high  energy  events  of  the  summer  of  1970,  phase-path 
fluctuations  data  were  collected  by  the  doppler  array  operated  by  Teledyne 
Isotopes.  Array  processing  of  the  data  was  accomplished  by  applying  co¬ 
phase  analyses  to  the  records.  Cophase,  a  statistic  characterizing  sig¬ 
nal  strength  as  a  function  of  velocity  and  direction,  has  been  defined 
by  Posmentier  and  Hermann  (1971). 

An  80-minute  time  window  was  moved  progressively  in  time,  in  steps 
of  15  minutes  starting  3  hours  before  and  ending  3  hours  after  the  time  at 
which  infrasonic  signals  were  expected  over  the  array.  The  frequency  band 
covered  by  the  analyses  was  from  0.000625  Hz  to  0.0020833  Hz,  corresponding 
to  periods  between  8  minutes  to  26.6  minutes.  The  width  of  the  frequency 
domain  window  was  0.0004166  Hz. 

A  sample  of  the  doppler  records  corresponding  to  some  of  the  events 
is  shown  in  Figure  22.  The  arrows  in  the  figure  correspond  to  the  arrival 
time  estimates.  Examination  of  this  figure  indicates  that  the  phase-path 
fluctuations  have  periods  between  20  minutes  and  60  minutes.  These  periods 
are  well  within  the  atmospheric  gravity  wave  period  range.  However,  the 
results  of  the  array  processing  showed  that  high  values  of  cophase  occur  at 
phase  velocities  and  azimuths  which  do  not  agree  with  those  that  would  be 
expected  from  the  location  of  the  high  energy  events,  i.e.,  450-650  msec  1 
and  238C115°  respectively.  These  results  are  not  difficult  to  understand 
since  the  estimated  arrival  time  for  the  long  period  atmospheric  gravity 
waves  is  about  2300  UT.  It  has  been  shown  in  Section  4  that  at  the  time 
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around  2300  UT  the  background  motion  increases  substantially.  Hus  is 
illustrated  in  Figure  23  where  it  can  be  seen  quantitatively  that  the  powe 
spectra  of  the  background  motions  is  an  order  of  magnitude  larger  at  this 
time  of  day  than  it  is  during  the  previous  hours. 

This  behavior  of  the  background  makes  the  identification  of  the 
signals  coming  from  high  energy  events  almost  impossible.  Ilte  direction 
of  arrivals  and  the  phase  velocities  obtained  from  the  analyses  suggest 
that  the  fluctuations  are  dominated  by  TID  activity. 


4*  J-m-nuJ.?  V*1*  J"  ****1^1  .  j;  |.,\  .*-  .  i'  r  .-  ||f;;  p 
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7.  SUMMARY 


Analyses  of  atmospheric  wave  data  from  a  long  period  multi-sensor 
airay  were  carried  out.  The  array  was  situated  in  the  New  York-New  Jersey 
region  and  operated  during  1969-1970.  The  results  of  the  analyses  enable  us 
to  draw  the  following  conclusions. 

Saturn  Apollo  launches  generate  acoustic  gravity  waves  in  the  fre¬ 
quency  range  0  00333  lit  -  0.0111  Hz.  These  waves  induce,  at  ionospheric 
heights,  phase-path  (doppler)  fluctuations  which  can  be  detected  by  phase 
sounders  at  distances  of  the  order  of  1S00  km  from  the  rocket  trajectory. 

The  direction  of  arrival  of  the  signals  intersects  the  trajectory  of 
the  rockets  at  a  point  at  which  the  altitude  of  the  rockets  is  approximately 
160  kn  Hits  fact  suggests  that  the  atmospheric  disturbances  are  generated 
at  ionospheric  levels. 

Seasonal  variations  m  the  lower  atmosphere  do  not  affect  these 
arrivals,  indicating  that  propagation  of  these  waves  takes  place  in  the  upper 
atmosphere 

The  phase  velocities  of  these  signals  are  or  the  order  of  720  msec* 
while  the  group  velocities  range  from  220  msec**  to  450  msec’1.  The  average 
horizontal  wavelengths  are  of  about  130  km. 

The  lack  of  arrivals  at  times  corresponding  to  smaller  rocket  launch¬ 
ings  and  the  wavelengths  involved  strongly  suggest  that  the  waves  are  generated 
by  rocket  engine  exhaust  plume  rather  than  by  a  ballistic  mechanism. 

Analyses  of  ground  level  pressure  data  have  tailed  to  detect  any  sig¬ 
nals  corresponding  to  the  doppler  signals.  This  lends  additional  support  to 
the  hypothesis  that  propagation  of  these  waves  takes  place  at  ionospheric 
heights. 
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A  dispersion  study  of  a  TIL)  was  used  to  illustrate  the  potential 
of  this  type  of  study  to  provide  information  about  the  neutral  stiucture 
of  tin:  upper  atmosphere.  The  horizontal  phase  velocities  of  this  T1D  event 
ranged  from  120  msec-1  to  240  msec"1  A  criterion  given  by  Tolstoy  (1971) 
was  used  to  test  the  proximity  of  the  observation  point  to  the  turning  point 
of  the  internal  gravity  wave,  and  then  a  WKB  approximation  was  used  to  com¬ 
pute  the  various  modes 

It  is  possible  to  construct  parametric  models  to  describe  the  time 
variation  of  the  ionospheric  background  motions.  It  was  found  that  in 
motions  the  frequency  range  1  to  13  CPH  respond  to  an  excitation  mechanism 
having  penodicties  of  12  and  24  hours.  Between  36%  and  50%  of  the  variance 
can  be  explained  by  the  model. 

A  simple  model  for  computing  ionospheric  response  to  internal  gravity 
waves  showed  that  for  horizontal  wavelengths  of  the  order  of  400  km,  at  a 
magnetic  dip  angle  of  72^,  the  amplitudes  of  the  doppler  variations  can  vary 
by  a  factor  of  3  depending  on  the  direction  of  the  wave  For  smaller  wave¬ 
lengths  this  anisotropy  is  less. 

Finally,  the  analyses  of  data  corresponding  to  high  energy  events 
showed  that  at  the  expected  arrival  time  (2300  UT)  the  doppler  backgiound 
level  increases  by  an  order  of  magnitude,  making  the  identification  of  the 
signals  almost  impossible. 
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